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Owing to the presence of the Coulomb barrier at astrophysically relevant kinetic energies it is 
very difficult, or sometimes impossible, to measure astrophysical reaction rates in the laboratory. 
That is why different indirect techniques are being used along with direct measurements. Here we 
O . address two important indirect techniques, the asymptotic normalization coefficient (ANC) and the 

D ■ Trojan Horse (TH) methods. We discuss the application of the ANC technique for calculation of 

^0 ' the astrophysical processes in the presence of subthreshold bound states, in particular, two different 

mechanisms are discussed: direct capture to the subthreshold state and capture to the low-lying 
bound states through the subthreshold state, which plays the role of the subthreshold resonance. The 
ANC technique can also be used to determine the interference sign of the resonant and nonresonant 
(direct) terms of the reaction amplitude. The TH method is unique indirect technique allowing 
one to measure astrophysical rearrangement reactions down to astrophysically relevant energies. 
' We explain why there is no Coulomb barrier in the sub-process amplitudes extracted from the TH 

£f~) , reaction. The expressions for the TH amplitude for direct and resonant cases are presented. 

o 

Q\ • PACS numbers: 26.20.+f, 21.10.Jx, 25. 55. Hp, 27.20.+n 

O ' 

in 

O ' I. INTRODUCTION 



For better understanding stellar evolution, cross sections of astrophysically relevant nuclear reactions should be 
known at the Gamow energy with an accuracy better than 10% Q. The presence of the Coulomb barrier for colliding 
^ ' charged nuclei makes nuclear reaction cross sections at astrophysical energies so small that their direct measurements 
in laboratories is very difficult, or even impossible. That is why direct measurements are being done at higher 
energies and then extrapolated down to the Gamow energy. Such an extrapolation procedure can cause an additional 
. i-H , uncertainty. Also for nuclear reactions studied in laboratory, the electron clouds surrounding the interacting nuclei 
lead to a screened cross section which is larger than the bare nucleus one (see Q,C|QLEi an( ^ references therein) . The 
enhancement factor is determined by the electron screening potential which is a model dependent quantity and its 
value in the laboratory is different from the one present in the stellar environment. There are four often used indirect 
techniques: the Asymptotic normalization coefficient (ANC) method 0, Coulomb breakup processes 0,111], Trojan 
Horse (TH) H0 and the Surrogate reactions method (see [l0( and references therein) . In this work we address only 
two indirect techniques, the ANC and TH methods. 



II. ANC METHOD 

The ANC method has been suggested in 0,^2 an d can be used to determine the astrophysical factors for peripheral 
radiative capture processes. The method can be applied for analysis of direct radiative capture processes leading to 
final loosely bound states. Due to small binding energies and strong Coulomb barrier, the direct capture reactions 
are peripheral. In previous papers jlll Il2l Il3| it has been pointed out that the overall normalization of the cross 
section for a direct radiative capture reaction at low binding energy is entirely defined by the ANC of the final bound 
state wave function into the two-body channel corresponding to the colliding particles. The ANC technique turns out 
to be very productive for analysis of the astrophysical processes in the presence of the subthreshold state . Here 
we address some applications of the ANC method in the presence of the subthreshold state. We also demonstrate 
how ANC technique can be used to determine the interference sign of the direct and resonant amplitudes for some 
important astrophysical radiative capture reactions. 
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A. Definition of the ANC 

We present first some useful equations for the ANC. Let us consider a virtual decay of nucleus c into two nuclei a 
and b. First we introduce the overlap function I of the bound state wave functions of particles c, a, and b |l5j : 

J ab(r) = < <Pa((a)<fb((b)\tPc((a,(b;r) > 

= i lc < JaM a j c m jc \J c M c > 

l c m lc j c m jc 

x < J b M b l c miJjcm jc > Yi cmic (r)Il Uejc {r). (1) 

Here (pi, Ci, Ji and Mi are the bound state wave function, a set of internal coordinates including spin-isospin variables, 
spin and spin projection for nucleus i. Also r is the relative coordinate of the centers of mass of nuclei a and b, 
r = r/r, j c , rrij a are the total angular momentum of particle b and its projection in the nucleus c = (ab), l c , mi a 
are the orbital angular momentum of the relative motion of particles a and b in the bound state c = (ab) and its 
projection, < ji77iij 2 m-2 1^3^-3 is a Clebsch-Gordan coefficient, Y"; cmc (f) is a spherical harmonic, and Iabi a j c ( r ) * s ^ ne 
radial overlap function which includes the antisymmetrization factor due to identical nucleons. The summation over 
l c and j c is carried out over the values allowed by angular momentum and parity conservation in the virtual process 
c — > a + b. The asymptotic normalization coefficient C c ahl • defining the amplitude of the tail of the radial overlap 
function I^ bl . (r) is given by [l^ 

jc i s r>R N c W_ v ^i c + 1/2 {2K ab r) 

1 abl a j c V) * L/ abl c j c ~ ' \ z ) 

where Rn is the nuclear interaction radius between a and b, W - Vc ,; c +i /2(^abf) is the Whittaker function describing 
the asymptotic behavior of the bound state wave function of two charged particles, k = \J1 [i a b £c is the wave number 
of the bound state c = (ab), fi a b is the reduced mass of particles a and 6, e c is the binding energy of the bound state 
(ab) and r\ a = Z a Z^e 2 Hab/K is the Coulomb parameter of the bound state (ab), Zi e is the charge of particle i. We 
use the system of units such that K = c = 1. There is another definition of the ANC, the most model independent 
one. The elastic a + b scattering amplitude in the channel (l c , j c ) has a pole in the momentum plane [T^ 

M lcjc (k) = SlcJc : - 1 h ^ p Wlcj ° . (3) 
2 1 k 2 % kip k kp 

corresponding to the bound state c = (ab) for k p = i k and to the resonance for k p = kn, where fcfj = fco — i ki is the 

resonance location in the momentum plane. Here, Si c j c is the elastic matrix element of the S-matrix. The residue in 
the pole Wi ajc is 

Wi cjc = ie^ (C c ablcj f, kp — i k, (4) 

Wi cjc =-(-l) h i(C c aUcMR} ) 2 , k p = k R . (5) 

For narrow resonances, ki « ko, 

(C c ablejc(R) ) 2 = rio . c . (6) 

Here r/Q is the Coulomb parameter for the resonance at momentum fco, ^z c j c (fco) is the potential (non-resonant) 
scattering phase shift taken at the momentum fco. Thus the residue in the bound state or resonance pole is expressed 
in terms of the ANC and for the resonance the ANC can be expressed in terms of the partial resonance width [l4| . 
Note that Eq. © holds only for k in the closest vicinity of the pole. For elastic scattering at positive energies in the 
presence of the Coulomb barrier, the elastic scattering amplitude with the bound state pole behaves (in the R matrix 
approach) as 

Mt , (fc) ^= -— e -*(*io-°i„) 1£ (7) 

JcV ' 2k E + e c + iT c /2 V ' 

where 



r c = 2P lc (E)j 2 c . 



(8) 
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Here Pi c (E) is the penetrability through the Coulomb-centrifugal barrier, 4>i c is the solid sphere scattering phase shift 
in the partial wave l c and cr/ e = tan _1 ( — ) , tq is the channel radius, 7^ is the effective (observable) reduced 

n=l 

width: 



7c = 2^ 7 (<W>) ■ (9) 



Thus at positive energies, E — > +0 due to the presence of the Coulomb-centrifugal barrier the elastic scattering 
amplitude behaves as the resonant scattering amplitude with the resonance width expressed in terms of the ANC. 
At positive energies the elastic scattering cross section in the presence of the bound state and the barrier behaves as 
the high-energy tail of the resonance located at energy E = — e c . That what is called the "subthreshold" resonance. 
However, it is not a resonance because the real resonance is located at complex energies on the second energy 
sheet, while the subthreshold resonance is just the bound state located on the first energy sheet at negative energy, 
corresponding to the bound state. At negative energies (positive imaginary momenta) Eq. reduces to Eq. 
Definitions of the ANC dictate the experimental methods of its determination. The ANC can be determined from 
peripheral transfer reactions which are dominated by the tail of the overlap function. Eq. © offers another possibility 
to determine the ANC, namely, by extrapolating the elastic scattering amplitude (or equivalently the phase shift) to 
the bound state pole jig. 



B. ANC and astrophysical processes 

(i) For peripheral direct radiative capture reaction a + b — > c + 7 to the final state l c j c proceeding through the EL 
transition, the cross section is 

° ~ \<i: b l cj M rL \Mr)>\ 2 



^jV ^' , ^ /a(W) lr £ |^, i (r)>l a . 



(10) 



Here L is the multipolarity of transition, tp^ i i (r) is the initial a + b scattering wave function with the relative 
momentum fc, in the partial wave U. Thus the ANC determines the overall normalization of the direct radiative 
capture cross sections. 

(ii) The elastic scattering amplitude Q describes the elastic scattering through the intermediate bound state c = (ab). 
Assume that it is an excited state. Then, when the excited bound state is formed it can decay into the ground state 
by emitting a photon. In this case we have the radiative capture process which is called the capture to the ground 
state through the subthreshold resonance. The amplitude of this process is given by 



I r l/2 r l/2 

M, 7 (k) -— e-^'c-*'.) c 7 -. (11) 



Here |ry 2 | 2 gives the radiative width for the transition from the excited bound state —> ground state. Thus in the 
presence of an excited bound state close to threshold, two different radiative capture processes can occur: direct 
capture to this excited bound state or capture to the low-lying bound states through this subthreshold bound state 
(capture through the subthreshold resonance) . In what follows we present some astrophysical reactions in the presence 
of the subthreshold state. 



C. ANC for 14 N + p -> ls O and the astrophysical S factor for 14 N(p,7) ls O 

The 14 N + p — > 15 + 7 reaction is a notorious example of an important astrophysical reaction where the sub- 
threshold state plays a dominant role. This reaction is one of the most important processes in the CNO cycle. As 
the slowest reaction in the cycle, it defines the rate of energy production jj and, hence, the lifetime of stars that 
are governed by hydrogen burning via CNO processing. The 14 N(p, 7) 15 reaction proceeds through direct capture 
to the subthreshold state 3/2 + , 6.79 MeV (binding energy 504 keV) and, possibly, via direct capture to the ground 
state and resonant capture through the first resonance and subthreshold resonance at E s = —504 keV. The over- 
all normalization of the direct capture is defined by the corresponding ANC. The ANC for the subthreshold state 
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FIG. 1: The 14 N( 3 He, d) 15 differential cross sections. The squares are data points and the solid lines are the DWBA calculations 
normalized to the experimental measurements in the main peaks; (a)- our data, (b) - our fit of the angular distribution measured 
in Ref. . 

E s = —504 keV also determines the partial proton width of the subthreshold resonance. In order to determine the 
ANCs for 14 N + p — > 15 0, the 14 N( 3 He, d) 15 proton transfer reaction has been measured at an incident energy of 
26.3 MeV. Angular distributions for proton transfer to the ground and five excited states were obtained. Angular 
distributions of deuterons from the 14 N( 3 He,<i) 15 reaction leading to the most important transition to the fourth 
excited state 3/2 + , 6.79 MeV in 15 measured by us at an incident energy of 26.3 MeV and in measured at an 
incident energy of 20 MeV, together with our DWBA fits are shown in Fig. ^ The proton ANC that we obtain for 
the 14 N + p -> 15 0(3/2+, 6.79 MeV) is C 2 = 27.1 ±6.8 fin -1 . Using our ANCs, we calculated the astrophysical factor 
and reaction rates for the 14 N(p, 7) 15 process. The capture to the 3/2 + , 6.79 MeV state dominates all others and 
the calculated astrophysical factor is 5(0) = 1.40 ± 0.20 keVb. The calculated and experimental S(E) factors for 
the transition to this subthreshold state are presented in Fig. [21 The uncertainty in 5(0) is entirely determined by 
the ANC of this state and the 13% systematic uncertainty in the experimental S(E) factor j5- We find that the 
astrophysical factor for the capture to the ground state is 5(0) = 0.15 ±0.07 keVb. The total calculated astrophysical 
factor at zero energy is 5(0) = 1.70 ± 0.22 keVb what is in excellent agreement with the 5 factor 5(0) = 1.70 ± 0.22 
keVb obtained from recent direct measurements performed at LUNA [l9j. The lower astrophysical factor of the 
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FIG. 2: The 14 N(p, 7) 15 astrophysical S factor for capture to the fourth excited state ((c): 3/2 + , 6.79 MeV), which includes 
the incoherent sum of the resonant and nonresonant terms. The squares are data points |l8|j: the solid lines represent the 
calculated S factor using our measured ANC. 

14 N(p, 7) 15 reaction leads to an increase in the age of the main-sequence turnoff in globular clusters pcj . 

D. ANC and interference of direct and resonant amplitudes 



To demonstrate how the information about the ANC can be used to determine the interference sign of the resonant 
and direct amplitudes of the radiative capture process we use the R matrix approach. Let us consider the radiative 
capture reaction a + b — > c + 7. 

The i?-matrix radiative capture amplitude to a state of nucleus c with a given spin J/ and relative orbital angular 
momentum of the bound state I / is given by the sum of resonant and nonresonant (direct capture) amplitudes |2l| : 



(12) 



Interference effects only occur in Eq. i|12|) if the resonant and nonresonant amplitudes have the same channel spin 
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/ and orbital angular momentum ij. In the one- level, one-channel approximation, the resonant amplitude for the 
capture into the resonance with energy E rn and spin J i; and subsequent decay into the bound state with the spin Jf, 
is given by 

" I J E-E rn + i^ ' (13) 

Here Ji is the total angular momentum of the colliding nuclei a and b in the initial state, J a and Jj, are the spins of 
nuclei a and b, and /, k, and li are their channel spin, relative momentum and orbital angular momentum in the initial 
state. Un i j f j i is the transition amplitude from the initial continuum state (Ji, I, li) to the final bound state (Jf, I). 
Also [TjliE)} 1 / 2 is real and its square, r j; \ (E), is the observable partial width of the resonance in the channel a + b 
with the given set of quantum numbers, [rA (E)] 1 / 2 is complex and its modulus square is the observable radiative 
width: 

r^^HF^)] 172 ! 2 - (14) 

The energy dependence of the partial and radiative widths is given by 

^u{E) = ^^T J ll (E Rn ) 1 (15) 

and 

^^)-(£^) 2i+lr ^(^«)' ( 16 ) 

respectively. Here, r^j .(-Erk) and T^j (Eji n ) are the experimental partial and radiative resonance widths, Sf is the 
proton binding energy of the bound state in nucleus A, L is the multipolarity of the gamma quanta emitted during 
the transition, and Tj i w y\ T^.. In a strict i?-matrix approach 

[^' Jf (E)} 1 / 2 ^2[P h (E)} 1/2 l J 7 j f - (17) 



Here the radiative reduced-width amplitude J^j f is given by the sum of the internal and external (or channel) 
reduced-width amplitudes: 

rjj f =itj f (mt) + ^ jf (ch). (is) 

Hence the total radiative width is 

\[^ Jf (E)]\ = \[T^ Jf (E)]lit + Kj+EtJuW (19) 

K*J, WfilU = 2 [ p h (E)] 1/2 7 7 % (int, ch) . (20) 
While the internal reduced- width amplitude is real, the channel reduced- width amplitude is complex [2l| and is defined 



as 



Z a e + ( _ 1)L Zte\ x , /(L + l)(2i+l) 



m L m% ) " V L 

(2L + 1)11 ^ a)L+V2 Cj ^ \l T ™S E ti 
([F h (k,a)} 2 + [G h (k,a)} 2 ) x W t f (2 k a) (l t L0\l f 0) 

U{Ll f J t I]kJf)JL{Ul } ). (21) 
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The nonresonant capture amplitude is given by 

TjNR _ _/o 1 3/2 1 l,+i-!,+l 1 Itl/2 

u IliJfJi — \ L ) 1 h " 



Z a e uL Z b e\ / (L+1)(2L + 1) 

X (2L+1)!! ( ^ a)L+V2 ro) 
xG u (k, r ) W_ Vftlf+1/2 {2Kr ) 

yJ~P]~{liO L0\lf0) U(L If Ji I; hJ f ) 

xj' L (hl f ), (22) 

^ {E) = Ff^J+Cf^roY (23) 

where i 7 i i and Gi t are the regular and singular (at the origin) solutions of the radial Schrodinger equation, k = ^J2ji^ej 

is the wave number, and fc 7 = E + ef is the momentum of the emitted photon. Integrals Jl{U If) and J L {h If) are 
expressed in terms of Fi i , Gi i and Whittaker function W 7 .^ ; /+1 / 2 and are given in 0,|23|- Both the channel radiative 
width and nonresonant amplitude are normalized in terms of the ANC, Cj f n f , which defines the amplitude of the tail 
of the bound state wave function of nucleus c projected onto the two-body channel a + b with the quantum numbers 
J/, /, If. Such a normalization is physically transparent: both quantities describe peripheral processes and, hence, 
contain the tail of the overlap function of the bound wave functions of c, a and b, whose normalization is given by 
the corresponding ANC. Note that in the i?-matrix method the internal nonresonant amplitude is included into the 
resonance term. Also, in the conventional i?-matrix approach the channel radiative width and nonresonant amplitude 
are normalized in terms of the reduced width amplitude, which is not directly observable and depends on the channel 
radius. However, it is more convenient to express the normalization of the nonresonant amplitude in terms of the 
ANC that can be measured independently [14j. Then only the radial matrix element depends on the channel radius. 
As we can see from Eqs (|21|l and <|22[l the relative phase of the channel radiative width and the nonresonant amplitude 
is fixed because only the ANC has unknown phase factor. Thus by measuring the ANC for the bound state we are 
able to fix the absolute normalization of the channel radiative width and nonresonant amplitude simultaneously. 



E. Interference of the resonant and nonresonant amplitudes for the 11 C(p, 7) 12 N astrophysical radiative 

capture reaction 

The evolution of very low-metallicity, massive stars depends critically on the amount of CNO nuclei that they 
produce. Alternative paths from the slow 3 a process to produce CNO seed nuclei could change their fate. The 
n C(p, 7) 12 N reaction is an important branch point in one such alternative path. At energies appropriate to stellar 
evolution of very low-metallicity, massive stars, nonresonant capture to the ground state and interference of the second 
resonance and the nonresonant terms determine the reaction rate. The ANC for 12 N — > n C +p has been determined 
from peripheral transfer reaction 14 N( n C, 12 N) 

13 C at 10 MeV/nucleon [23. The contributions from the second resonance and interference effects were estimated using 
the R-matrix approach with the measured asymptotic normalization coefficient and the latest value for the radiative 
width of the second resonance 23j. The ANC gives useful information not only about the overall normalization 
of the direct capture amplitude, but also about the radiative width of the resonances. According to Eqs. I|2(J|I . 
the channel part of the radiative width may be determined from the ANC. Since the channel part is complex, 

\P-yJf Ji {Ey\Jh = A + it, while the internal part of the radiative width amplitude is real, i. e. [T 7 j, j i (E)]l^ = v, the 
total radiative width is given by 

T lJfJi {E) = (\ + v) 2 +t 2 . (24) 

The relative phase of A and v is, a priori, unknown, so these real parts may interfere either constructively or destruc- 
tively. Hence, r 2 always provides a lower limit for the radiative width and additional stronger limits may be obtained 
if assumptions are made about the interference between the two real contributions. For constructive interference of the 
real parts, the channel contribution gives a stronger lower limit. In the case of destructive interference, if |A| > \v\, the 
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channel contribution gives an upper limit for the radiative width. These limits depend on only one model parameter, 
the channel radius. 

Recently, a measurement at RIKEN j23| found the gamma width to be 13.0 ± 0.5 meV. Using the measured ANC 
we find that for a channel radius of ro — 5.0 fm, r 7 j. j i (En) c h — 54 meV. Taking into account the experimental value 
of the total radiative width, one can find the internal contribution from 

T lJf j i (E R ) = \T lJf Jt (Er)]/ 2 + T lJf j t (E R )H 2 t \ 2 . (25) 

There are two solutions, 15 and 112 meV. Assuming that the second value is too high [24j . we conclude that the internal 
part of the radiative width is 15 meV, and destructive interference between the real parts of the channel and internal 
contributions gives the experimental value, 13 meV. In this case, the channel contribution alone represents an upper 
limit for the radiative width, while the square of the imaginary part of the channel contribution, 1.8 meV, gives a lower 
limit. The relative phase between the direct capture amplitude and the channel contribution to the radiative width 
of the second resonance is fixed in the R-matrix approach. Therefore, when the channel contribution to the radiative 
width dominates, the sign of the interference effects may be determined unambiguously. For n C(p, 7) 12 N, we find 
that the nonresonant and resonant capture amplitudes interfere constructively below the resonance and destructively 
above it. It has important consequences on the reaction rates for 12 N production. In particular, the reaction sequence 
7 Be(a, 7) 11 C(p, 7) 12 N 7 Be(a, j) 11 C(p, 7) 12 N will provide a means to produce CNO nuclei, while bypassing the 3 a 
reaction, in lower-density environments than previously anticipated |25| . 



F. Interference of the resonant and nonresonant amplitudes for the 13 N(p, 7) 14 astrophysical radiative 

capture 

13 N(p, 7) 14 is one of the key reactions which trigger the onset of the hot CNO cycle. This transition occurs 
when the proton capture rate on 13 N is faster, due to increasing stellar temperature (> 10 8 K), than the 13 N j3- 
decay rate. The rate of this reaction is dominated by the resonant capture to the ground state of 14 through 
the first excited state of (Er = 0.528 MeV). However, through constructive interference, direct capture below the 
resonance makes a non-negligible contribution to the reaction rate. We have determined this direct contribution by 
measuring the asymptotic normalization coefficient for 13 N + p — > 14 O(0.0 MeV). This ANC has been determined 
from the peripheral reaction 14 N( 13 N, 14 0) 13 C [2(|. The radiative capture cross section was estimated using an R- 
matrix approach with the measured asymptotic normalization coefficient and the latest resonance parameters. What 
is not known is the sign of the interference term between the resonant and nonresonant components of the radiative 
capture amplitudes. As we have mentioned it is possible to sometimes infer the sign of the interference to be used 
in an R-matrix calculations of the radiative capture cross section if the ANC is known even in the absence of direct 
experimental data. Such is the case for the reaction being considered here. At energies below the resonance, the 
channel part, which depends on the ANC, has the same sign as the nonresonant amplitude leading to the constructive 

1/2 

interference of these two terms. From Eqs. (|2*T|) and (|2*U|) we find \Pl f i j f (Er)] f = 0.90 + i 0.02 eV 1 / 2 and the channel 
radiative width \[T^ Jf (E R )} | = 0.81 x 10~ 6 eV at the resonance energy and the channel radius ro = 5 fm. The total 

resonance radiative width is \[T^' Jf (E)}\\ = \[T^ Jf (E)]l^ + [T^ Jf (E)] / \ . Thus there are two possible solutions for 
the internal part, a large negative value T (E)] ^ — —2.73 eV 1 / 2 and a small positive value IT'I 1 , (E)] ^ =0.93 

eV 1 / 2 . The first solution leads to the destructive interference with the non-resonant component at energies below 
the resonance, but it yields a high internal radiative width, \^ i j f (E)] < J = 7.48 eV. The second solution leads to 
the constructive interference with the non-resonant component at energies below the resonance peak. We select this 
second solution because it is corroborated by the microscopic calculations |22|, where it has been shown that the 
internal and external parts of the El matrix elements have the same sign and very close magnitudes. Our choice is 
also supported by the single-particle calculations j2|| [2^ . Due to this constructive interference we find the S factor 
for 13 N(p, 7) 14 to be larger than previous estimates. Consequently, the transition from the cold to hot CNO cycle 
for novae would be controlled by the slowest proton capture reaction 14 N(p, 7) 15 0. 



III. TROJAN HORSE 



The Trojan Horse method (THM) is a powerful indirect method which selects the quasi- free (QF) contribution 
of an appropriate three-body reaction performed at energies well above the Coulomb barrier to extract a charged 
particle two-body cross section at astrophysical energies free of Coulomb suppression. The THM has been suggested 



9 



by Baur [9j and has been advanced and practically applied by a group from the Universita di Catania working at 
the INFN-Laboratori Nazionali del Sud in Catania in collaboration with other Institutions (see [j| and references 
therein). The THM has already been applied many times to reactions connected with fundamental astrophysical 
problems [H such as 7 Li(p, a) 4 He, 6 Li(d, a) 4 He, 
6 Li(p, a) 3 He, and many others, see [B| and references 
therein. 

Let us consider the TH reaction 

a + y + b + B, (26) 

where a = (xy). The subreaction of interest is 

x + A^b + B. (27) 

In the TH method the incident particle a is accelerated to energies above the Coulomb barrier. After penetration 
through the barrier the projectile breaks into x + y leaving the fragment x to interact with target A, while the second 
fragment-spectator y leaves carrying away the excess energy. By a proper choice of the final particle kinematics, 
the THM allows one to extract the cross section of the sub-process (|27f) . However the extracted amplitude of the 
reaction (|2*7Jl in the THM is half-off-energy shell because the initial particle x in the sub-process <(T7|) is off-the-energy 
shell. It has been suggested in the original paper that the virtually of particle x is compensated for by the 
higher momentum components in the Fermi motion of the fragments x and y inside the projectile a. However, high 
momentum components means that the distance between the fragments is so small that the interaction between the 
fragments is not negligible and the mechanism of the reaction is more complicated than the QF one. Instead, the 
virtuality of particle x in the extracted cross section is significantly compensated if we take into account the binding 
energy of the fragments x and y in the projectile a |3lj . 

The THM allows one to determine both direct and resonant reactions (|27|l . As an example of the result achieved 
using the THM, we present in Fig. [31 trie astrophysical factor for the 3 He(<i, p) 4 He process determined from the 
3 He( 6 Li, ap) 4 He TH reaction [3^. The TH resonant cross section (full dots) is normalized to the direct experimental 
data (open circles and open triangles) at energies near the resonance peak. The black solid line is the result of a fit 
of the TH data (see ref.[33 for details), showing the trend of the bare nucleus S(E)-factor, while the blue solid line is 
obtained by interpolating the screened direct data. 



A. TH reaction amplitude 

A simple mechanism describing the TH process is the so-called QF process shown in Fig. 0] In the quasi-free 
process it is assumed that the incident particle (assume incident particle is A) interacts with one of the fragments of 
a = (xy), say with x which is considered to be "quasifree", while the second fragment is considered to be a "passive" 
spectator which is not involved in the process. Thus the interaction of the spectator y with x and A in the knockout 
process is neglected . The fact that the fragment x is not free is taken into account by folding the quasi-free reaction 
amplitude with the Fourier component of the (xy) bound-state wave function which takes into account the Fermi 
motion of x in the bound state a = (xy). 

In this section we present a derivation of the TH reaction amplitude from the general 2 — > 3 reaction amplitude for 
the TH process (I26[) . A general expression for the amplitude of the reaction is given by 

M =< x&xS <p v <Pb<p B \AVf(l + G+AVA^AcpaX^ > 

(28) 

=< X^xi'JwPBK^Vf G+ + l)AV i \<p A <Pax'? ) > ■ 

(29) 

The amplitudes (|28|l and l|29l) are the post and prior forms of the exact amplitude. Let us consider the post form. 
Here, G + is the total Green function of the system a + A, xi + ^ is the distorted wave describing the scattering wave 
function of a + A in the initial state of the reaction, \\,b 1& * ne distorted wave describing the scattering of particles 
b + B in the final state: the distorted wave XyF describes the distorted wave of the spectator y and the center of 
mass of the system F = b + B in the final state. For the moment we assume that Coulomb interactions are screened. 
Eventually we can take the limit of the screening radius to infinity. Also fa is the bound state wave function of nucleus 



m = V aA - U aA , 

AVf = V bB - U bB + Vy F - Uy F , 



(30) 
(31) 
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20 




E c . m . (MeV) 

FIG. 3: The 3 He(d,p) 4 He astrophysical S factor determined from the TH reaction. The open circles and open triangles are 
direct experimental data; the full dots are the TH data. The black solid line represents the behavior of the bare nucleus 
S(E)-factor, resulting from a fit on the TH data, while the solid blue line is interpolation of the direct data. 



Vij and Uij are the interaction potential and optical potential between particles i and j. For example, V a A = V x a+V v a- 
To extract the amplitude of the subprocess x + A — > b + B, which is the final goal of the TH method, we note that 
the Hamiltonian of the system a + A is 

H = H a A + H a + Ha = H x a + H y p + H x + Ha + H y , (32) 

where Hi is the internal Hamiltonian of nucleus i and fly = + Vij is the Hamiltonian of the relative motion of 
nuclei i and j, is their relative kinetic energy operator and Vij is their interaction potential. The total Green's 
function operator can be written as 

E-H aA -H a -H A +iO K ' 

1 



E-H xA -H yF -H° xyA + iO 



E - H xA - T y p - U yF - AV yF - H Q xyA + i0 



(34) 
(35) 
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a 



A 




FIG. 4: Pole diagram describing the quasi-free mechanism. 



G+ + G + AV yF G 



Here AVf = V yF - U yF , V yF = V yx + V yA , 
H °yA = H x + H v + H A and 

G+ = 



(36) 



(37) 



E-H xA - T yF - U vF - H° xyA + tO 
We substitute Eq. I|37|l into H28[) and drop the term 

AVf G + AVyF G + AVi as the higher order term in the perturbation expansion over AV. Then we get from Eq. 1281) 



M =< xUxyF^yVbVslAVfil + G+ AV^a^xT' > 



(+) 



(38) 



To single out the TH subprocess amplitude we replace AVi = V X A+V y A — U a A by V x a and AVf — Vyp — UyF+VbB — UbB 
by AVbB = VbB — UbB- Then the amplitude l|28ll becomes 



M =< x[ B ] XyF L Py l Pb l ~PB\AV bB (1 + G+ V xA )\ipA<Pa,X 



(+) 



> 



= < ^(-)J-) 



XbB XyF l Py l Pb<PB\AV b B (1 + G+ A V xA )\^>A^aXi 



(+) 



> 



Here 



G , 



1 



E x a - H xA + iO 



(39) 
(40) 



12 



X b 




interactions. 

and E x a is the relative kinetic energy of particles x and A. The appearance of G^ A in Eq. I|40|) is due to 

< X i y~F<*Py l PbtPB\G + =< X [ ~FVy<£btPB\G+ A . (41) 

Eq. (|39J) reveals a very important result. It contains a factor 1 + Gt A V X A- F° r the on-shell case, the relative 
momentum of particles x and A p x A — k x A, where k x A is the x — A on-shell relative momentum related with their 
relative kinetic energy as E x a = V x aI '(2 Hxa)- Correspondingly 

(1 + G + V xA )W^* A >= xt A (r xA ). (42) 

is the scattering wave function of particles x and A interacting via the optical potential V x a- We assume at the 
moment that all the Coulomb interactions are screened. However, in the TH method the entry particle x is not free 
because it is in the bound state a = (xy), i. e. the momentum of x is not fixed. In other words, x is off-the-energy 
shell because E x a ^ p xA /{2 ^ x a)- For the off-shell case 

(l + G+^)|e*P— >=xt os) u xA , P ^A) (43) 

is the so-called off-shell scattering function, 

B. TH method for direct reactions 

We first consider the direct subreaction (|27|l . We assume that this reaction proceeds through the transfer of particle 
c from A to x (it can be also considered as a particle transfer from x to A), i. e. A = (Be) and b = (xc). The "pole" 
diagram corresponding to the on-shell reaction describing the particle c transfer mechanism with the x — A rescattering 
in the initial and b—B rescattering in the final state is shown in Fig. This diagram describes the DWBA amplitude. 
To simplify Eq. (|39|l in the case of the direct transfer subprocess, we insert the projection operators ^2 fx >< fx, 
^2<Pb >< fB and ^2 >< </?c into the bra and ket states. The sum is taken over discrete states and an integral is 
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used for the continuum states of the corresponding nucleus. We leave in the projection operator only the ground state 
projections tp x >< ip x , (p B >< <£b and ip c >< ip c assuming that only the ground states of x, B and c contribute to 
the reaction. If necessary the excited states can also be taken into account. Then we get 

< tp y (p b (p B \AV bB (1 + G+ A V xA )\(p A tp a > 
«< ipb\lfcfx >< (fi B \AV bB \ip B > 

x(l+ < ip x \G+ A \tp x >< <p x \V xA \<Px >) 

< VcfB^A >< VxPyl^a > (44) 

We introduce the overlap functions JSL =< tpp ifj\ip a > and use the approximation < (px\V xA \(p x > ~ U xA \ also we 
use the approximation 



< <p x \G+ A \<p x >« G x u } + = (E xA - T xA - U xA + i0)~\ 



(45) 



Note that < (p B \AVbB\(pB >~ V x b + V c b ~ UbB, where Vjb is the interaction potential between the point like nuclei 
j = x,c and B. All the neglected terms are higher order in the perturbation theory over AV. Then we get in lowest 
order for the TH amplitude with the subprocess described by the direct transfer reaction i|27|) : 

M =< X$ \x$%c\Wb\ IcB (1 + G + xA V xA )\ I a xyX { + ] > • 

(46) 

The expression in the brackets is the amplitude of subreaction (|27|l which is the final goal of the TH. To see it we just 
rewrite (14611 in momentum space: 

M = l^^4F-\PvF)M-% b B,P xA )I^ y (p xv ) 

xx! +) (p,a), (47) 

where 

m yPx - m x p y m y 

Pxy — j — Pa - Py (48) 

m x + m y m x 

Also note that in the center of mass of TH reaction a + A^y + b + B the relative momentum is given by p aA = p a 
and p y F = Py We denote by Pi (kj) the momentum of the virtual (real) particle i and by (ky) the relative 

momentum of virtual (real) particles i and j. Also XiiPxA) = X^\{Pxa), i- e. it is the Fourier component of the 
a — A scattering wave function with the incident momentum \i. a A which in the center of mass of the TH reaction is 
just k a . Correspondingly X y ~ F (PyF ) = Xk~t(Pyf )• 

The half-off-the-energy shell amplitude of the subprocess l|27|) is given by 

M™\k bB ,p xA ) =< X&ll\*V bB \I&xU^, PaA > ■ 

(49) 

The virtuality of the entry particle x of this amplitude results in the fact that the relative momentum of particles 
x and A in the initial state of reaction l|27|) p xA ^ \/2 [i xA E xA . Due to the off-shell entry particles amplitude l|49|) 
does not have the Gamow penetration factor. We would like to underscore that from E xA + Q = E b s for positive 
Q > for reaction (|2~7Jl at E xA — > 0, E b s ~ const. Hence the off-shell scattering function xf 08 )k A P A IS * ne 0Iu y ^ xA 
dependent factor in M sub at E xA — > 0. The off-shell scattering function is a universal factor which does not depend 
on the specifics of the direct reaction. Rewritting matrix element in Eq. (|49|l in the momentum representation gives 



M (h bB ,PxA) = J J2^ x ^b W bB > 

xl*j(p' x - ^ p 6 ) AV bB I* B {VB - ^ P A ) 



m b m A 
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Approximation AV b B ~ V c b, which works for m x > m c , is enough for us to investigate the dependence of 
M sub (k b B, Pxa) on E x a for arbitrary masses of x and c. Using this approximation we get from Eq. I|5U|I 



M sub (lc r. ^ - f dPbB dPxA (r> \ 

i k bB,PxA) - J J^X kbB (PbBj 

xKdPx Pb Wf B (p B P A ) 

m b m A 

X xtos)^ A , PscA (PxA)- (51) 

Here W^ B (p c s) is the form factor determined by 

W c A b (Pcb) = J dr^e-'P— U cB (r cB ) 7 c A B (r cB ). (52) 

The Fourier component of the off-shell scattering function X(o S )k A p a { v xa) is given by 

X( os )^ AiVxA {Pxa) = Kp'xA-Pxa)+G+{p xA \E x a) 

xT{p' xA ,Pxa;E xA ), (53) 

G+(p' xA ; E xA ) = — -, (54) 

E X A -p xA l^^xA + «0 

T(p x Ai PxA', E x a) is the off-shell x — A scattering amplitude. Amplitude M sub (k b B, Pxa) extracted from the THM 
should be compared with the corresponding on-shell reaction amplitude 

shn i \ f dPbs dp x A *(-), \ 
{k bB ,k xA ) - J J^E^k bB ^ bB > 

X-C KPx Pb AVbs Icb(pb Pa) 

m b - m A 

*X^a{Vxa) (55) 

Eqs. (|47(l and (|49|) is our final result. The diagram corresponding to this amplitude l|47|) is shown in Fig. Eq. I|47|) 
is a general expression for the TH reaction amplitude which contains the half-off-shell direct subprocess amplitude 
and the initial and final state rescatterings. As we can see the subprocess amplitude is not factorized, but instead is 
folded with the initial and final state distorted waves and the overlap function for a — ► y + x. Note that if the initial 
and distorted waves in the momentum space are replaced by delta-functions, Eq. (|4T|) just becomes a trivial plane 
wave impulse approximation described by the diagram of Fig. [3] 



M c 



C. TH for resonant reactions 



In Subsection IIII Al we derived a general expression, Eq. I|39l) . for the amplitude of the TH reaction l|2t)|) which is 
valid for both direct and resonant subprocesses J2ZJ). Here we consider the resonant TH reactions, i. e. we assume 
that the subprocess (|27|l proceeds through the intermediate resonance F* . Our goal is to relate the half-off-shell and 
on-shell resonant amplitudes. Note that it is easier to relate the off-shell and on-shell resonant reactions than the 
direct ones. The resonant TH amplitude can be extracted from Eq. I|39|) in a straightforward manner because it 
contains the Green's operator G^ A - Below we demonstrate how to do it. For simplicity here we neglect the initial 
and final state interactions. 

The TH amplitude of the reaction (|26[) . which proceeds through the resonance state F* in the intermediate system 
x + A, is given by 

M = M sub ^(k bB ,p xA ) I xy {Pxy)- (56) 

Here M sub ^ is the amplitude of the resonant subprocess (|27|) . Usually in practical calculations the overlap function 
I xy is expressed in terms of the corresponding single-particle bound state wave function ip xy : 

I xy = S^y 2 (p xy . (57) 
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FIG. 6: Diagram describing the TH reaction a + A^y + b + B proceeding through the direct subprocess x + A — > b + B 
mechanism. Bubbles show initial and final state interactions and the off-shell scattering function. 



Here, S xy is the spectroscopic factor of the bound state (xy) in a with given quantum numbers. For simplicity we 
don't write down symbols corresponding to the quantum numbers and assume that S xy = 1. In the momentum space 
the bound state wave function is given by 



/ W(Pxy) 
tPxyyPxy) — ^^xy 2 , J 



(58) 



W(Pxy) 



dre lp *v r V xy {r) y xy (r) 

,2 



Now we can find the virtuality factor 



(-£a 



E x 



E 



xy 



2/i 



! (PxyiPxy), 



Irrir, 



(59) 
(60) 

(61) 



of the virtual particle x using the energy and momentum conservation laws in both vertices a — > x + y and x + A — ► F* 
After simple algebraic transformations we get 



&x E X A 



PxA 



1 



2(1. 



xA 



■[Ply + «y) 2 } <0. 



(62) 



Thus we derived a very important result for the relative kinetic energy of particles x and A E x a in the TH method: 
E x a < Pxa/^^xA, i- e. always k x A < PxA, where k x A = \/2 [i x a E x a is the x — A relative on-shell momentum. The 
half-off-shell resonant reaction amplitude in the TH method is described by the diagram shown in Fig. [7] and is given 
by 



M s ^ R \k bB , PxA ;E) = ~i(47r) 2 J—^ 

2 V MbBKfcs 
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a 



A 




FIG. 7: Diagram describing the resonant reaction a + A^y + b + B. 



Here k. 



mo=— (o 



V'rfcB (-EbB , r ) Wi (p xA , k xA (R)) 



E x a - E [ X R 2 



MR) ~ V ^ l lxA ^xA ' kfcB is the on-shell relative momentum of particles b and _B in the final state, /o( m o) is 
the resonance orbital angular momentum (its projection), Yi omo is the corresponding spherical harmonics, r = r/r, 
Sfi g is the nonresonant (potential) scattering phase shift of particles b and B in the final state. The off-shell form 
factor 



(63) 



m {PxA,k xA{R) ) = I drr 2 ^(r)V(r)ji {p xA r) 



= (E l f2-E PxA ) / drr 2 ^(r)j lo (p xA r) 



{E x R 2-E PxA )^l{ PxA ) 



(64) 



Here ip^(r) is the resonant Gamow radial wave function, ^ni f '(p xA ) is its Fourier component, ji (p xA r) is the 
spherical Bessel function, E PxA = P xA /2 ^ xA , n is the principal quantum number. 



(R)i 
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Let us write down the well known expression for the on-shell Breit-Wigner resonance amplitude for the resonant 
process x + A — ► b + B 



M^(k bB , k xA ;E) = --Unf J— } 



1 



4 V VbBhB V VxAkxA 

m =l 



y/T bB (E bB ,r ) y/T bB {E x A,ro) 

x p 



(65) 



where k^ is the on-shell relative momentum of the initial particles x and A and k^ is the on-shell relative momentum 
of the final particles b and B. In the i?-matrix method the resonance width contains the Coulomb-centrifugal barrier 
penetrability factor which exponentially decreases with energy. Hence for E x a — * the resonant amplitude M R ~ 
y/P^JkxA) M . Just this factor makes it difficult or impossible to measure resonant reactions at astrophysically 
relevant energies. Now we compare the half-off-shell resonant amplitude, Eq. i|63[l . aud the on-shell amplitude, Eq. 
(|65|1 . The half-off-shell amplitude contains the form factor 

wi (pxA, k( x A)R)- The barrier factor should come from the integral representation in Eq. I|64|) . namely from ji (j> x A r )- 
However, ji (jp x at) does not contain the Coulomb penetration factor and does not depend on the on-shell momentum 
k x A- Hence in limit k x A — > the off-shell form factor does not go to zero. We underscore that it is very important 
that always in the TH reaction p x A > k x A- Comparing Eqs 163|) and (|65[l we get 



M( s »(kiB, k xA ;E xA ) = -\e l5 W^ A \[—^- 

2 V fJ-xAk 



xA 

x ^Y [E l A,rQ \ MSub{R) ^ Pxa;E xA ). (66) 

Wl {p x A,k xA (R)) 

Note the only difference between the half-off-shell and the on-shell resonant amplitudes is the appearance of the form 
factor wi (p x A, k X A(R) ) ■ Now we give the expression for the on-shell resonant cross section which can be derived from 
the TH half-off-shell resonant cross section 

frl \ HxAk x A (J-hBkhB 1 1 f ,_, 

<Exa) = — <m* — Wa^I kbB 



dn kxA \M iR Hk bB , k xA ;E xA ) | 2 (67) 

1 VbBhB 1 1 T xA (E X A,r ) 
4 (2tt) 2 k 2 xA 4n\w lo {p x A,k xA (R))\ 2 



y"dfi kbB J &n^ A \M suh( - R HUB, p x a;E xA 



)l- (68) 



IV. SUMMARY 



In this work we have addressed two important indirect techniques in nuclear astrophysics use, asymptotic normal- 
ization coefficient (ANC) and Trojan Horse (TH) 

method. Both techniques allow one to determine the astrophysical factors at Gamow peak or even at zero energy 
avoiding extrapolation procedure. The ANC method determines the overall normalization of the peripheral radiative 
capture processes. The ANC technique becomes especially powerful for astrophysical processes proceeding through 
a subthreshold state - a loosely bound state. In this case the ANC determines both the overall normalization of the 
direct radiative capture to the subthreshold state and the resonance partial width for captures through the subthresh- 
old resonance. We demonstrated the application of the ANC technique for the key CNO cycle reaction 14 N(p, 7) 15 0. 
The ANC method turns out to be useful also for determination of the sign of the interference term of the resonant and 
nonresonant radiative capture amplitudes. We demonstrated it for two important CNO cycle reactions: 11 C(p, 7) 12 N 
and 13 N(p, 7 ) 14 0. 

The TH method allows one to determine the astrophysical factors for astrophysical reactions, both direct and 
resonant. In practical applications the astrophysical factor extracted from the TH reaction is available in a wide 
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energy range from astrophysical energies to higher energies. Its absolute normalization is determined by normalization 
of the TH astrophysical factor to the one obtained from direct measurements at higher energies. Assuming that the 
energy dependence of the TH astrophysical factor is correct, one can determine the absolute astrophysical factor at 
astrophysical energies. In this work we have derived a general expression for the TH reaction amplitude which takes 
into account the off-shell effects and initial and final state interactions. The direct and resonant TH reactions are 
considered separately. We derived the TH amplitude for direct subreactions in terms of the off-shell scattering wave 
function. The energy dependence of this wave function determines the energy dependence of the TH astrophysical 
factor for an arbitrary direct reaction mechanism. We connect the TH resonant cross section with the on-shell resonant 
cross section. We intend to use the derived equations to calculate the absolute astrophysical factors. 
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